Vibration-induced granular segregation: a phenomenon driven by three mechanisms 
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The segregation of large spheres in a granular bed under vertical vibrations is studied. In our 
experiments we systematically measure rise times as a function of density, diameter and depth; for 
two different sinusoidal excitations. The measurements reveal that: at low frequencies, inertia and 
convection are the only mechanisms behind segregation. Inertia (convection) dominates when the 
relative density is greater (less) than one. At high frequencies, where convection is suppressed, 
fluidization of the granular bed causes either buoyancy or sinkage and segregation occurs. 

PACS numbers: 46.10.+Z, 64.75.+g, 83.80.Fg 
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Many theoretical and experimental studies have been 
carried out in the last five decades aimed to reveal 
the physics of one of the most intriguing phenom- 
ena in granular matter: vibration-induced segregation 

QBiaiiaiiliailllllllllll. However, 
when this phenomenon intensively investigated appeared 
to be well understood, new scientific puzzles came into 
scene [TKllllllllliaiSllllil. Air-driven segrega- 
tion, inertia, condensation, are now new words added 
to the already vast list of concepts important in the 
subject. Thus, since this problem is an important con- 
cern to industries dealing with granulates, these recently 
disclosed effects should be further investigated. Gran- 
ular segregation was first reported in 1939 by Brown 
|23l | and studied ever since by the engineering commu- 
nity Q, II S S 0, until it was brought in 1987 to the 
physics realm with the suggestive name of "Brazil Nut 
Problem" (BNP) Results related to this problem 

established themselves as benchmarks of granular segre- 
gation. But the question: why the Brazil nuts are on 
top, seems to be yet an open matter of discussion. Both 
theoretical and experimental studies have focused on the 
infiuence of size, friction, density and excitation param- 
eters II i, i, lla. Ill 111 El IH 111 113 and the results 
explain, or obscure, bit by bit the underlying mechanism 
behind the BNP. Some of these results support the idea 
that it is "void-filling" beneath large ascending particles, 
the mechanism promoting the upward movement of an in- 
truder in a shaken granular bed 14j . Other researches 
claim that global convection is the driven force behind the 
BNP ig], and others, that arching ^ or inertia [H 112 
are crucial elements to explain it . The dilemma is not 
yet settled on the verge of even more recent findings ; 
being the most relevant the surprising result that de- 
creasing the density of the intruder does not necessarily 
mean a monotonic increasing of the rise time, as might 
be previously suggested by studies in 3D and 2D [ig . 
Furthermore, based on computer simulations Hong et al 
even dared to predict the reverse segregation effect in the 
BNP iliiiial (known now in the literature as the RBNP), 
which was immediately confronted by two groups |20ll2^ . 
but nevertheless observed in the laboratory by Breu et al 
[25| . Finally, Yan et al (21(| recently failed to confirm the 



experimental findings of Mobius et al ^3 ■ 

Based on this debate, a simple, yet overwhelming con- 
clusion arises: more research is needed if we want to un- 
cover the physics of this elusive granular matter problem. 
This Letter aims to contribute to its final understanding. 

Our experimental set up consists of a Plexiglas cylinder 
(closed at one end) of 10 cm inner diameter and 26 cm 
of length. The Plexiglas cylinder is fixed to a vibrating 
table fed with an amplified periodic voltage coming from 
a function generator (HP-33120A). In a typical experi- 
ment, the column is filled with small seeds or glass beads. 
Rise times are measured by a stopwatch. Excellent re- 
producibility of the data is obtained if the temperature 
and humidity do not change during the experiments. 

In Fig^ we show the rise time as a function of the 
intruder relative density {pr = p/ Pb, where pb is the den- 
sity of the bed particles) for 7 different starting depths: 
5,7,9,11,13,15 and 17 cm below the surface, at 5 Hz and 
r = 3. The bed column had 21 cm of height (hence, at 17 
cm of depth our 4 cm diameter intruders touch the bot- 
tom of the container). The size and density of the bed 
particles we used (tapioca monodisperse spheres) are re- 
spectively 3.1 mm and 0.57 g/cc. The intruders are plas- 
tic spheres filled with different materials to change their 
densities. 

Our first clear observation is that the ascension dy- 
namics of the spheres whose starting positions were at 
the bottom of the column, has a monotonic dependence 
on Pr] the curve diverging at ~ 1- At any other 
depth, the spheres regardless their density segregate to 
the surface following a non-monotonic dynamics. This 
non-monotonic ascension dynamics was previously ob- 
served by the group of Chicago although the peak 
they observed was positioned at a relative density less 
than one and their measurements correspond to only one 
depth (around 5 cm). We can normalize the data with 
> 1 in Fig^by making the rise time of the heaviest 
sphere equal to one (for each depth). In this way, the 
rise time of every sphere will be measured with the time 
scale of the fastest (the densest) one. These normalized 
results, for the spheres denser than the granulate, col- 
lapse into the same curve, see Fig|21 Based on this result 
we believe that these spheres segregate mostly by inertia. 



FIG. 1: Dimensionless rise times for intruders with pr > 1 for 
all depths, normalized with the rise times of the fastest (dens- 
est) sphere, at each depth. The line is the best fit obtained 
with apr, giving a = 9.84 and b = —1.09. Inset: dimension- 
less rise times for intruders with pr < 1 for different depths, 
normalized with the rise time of the fastest (lightest) sphere, 
at each depth. 



This concept has been evoked in the hterature for some 
time already 0, El El El El , but never used to quan- 
tify granular segregation in 3D. In a recent paper |22j . 
we proposed a simple theoretical model to explain the 
rise dynamics of heavy spheres in a vibrated granulate. 
The model is based on energy considerations and states 
that, on each cycle, the kinetic energy of the intruder 
is lost by friction during its penetration into the granu- 
lar bed; l/2mVfg = (3{h)Pi (where vto is the "take-off 
velocity the bead has when the granulate reaches a neg- 
ative acceleration a = —g, m the mass of the bead, and 
P{h) the friction force exerted upon it by the granulate). 
Therefore, since in our experiments the volume of the in- 
truders and the vibration conditions are constants, the 
penetration length per vibration cycle should be directly 
proportional to the density. Thus, the number of cycles 
for a sphere to segregate to the surface is inversely pro- 
portional to it. The smooth line in Fig|21was obtained 
by using ap^, where the best fit gives b — —1.09. What 
is very interesting is that all peaks in Fig^are at ~ 1 
and afterwards, lighter intruders start to ascend faster. 
In the inset of Fig|2] we show normalized rise times for 
Pr < 1- For this case, normalization was done using also 
the fastest sphere (the lightest). We will come back to 
this plot later. 

We carried out a second experiment. The rise time of 
light spheres was measured again as we change pr, but 
now the spheres are always positioned at the bottom of 
a container filled by a granulate of 6 cm of height. The 
frequency we use is 50 Hz (with T = 3). We can see 
that the lighter the spheres the faster they ascend, see 
FigOl This is contrary to what happens at low frequen- 
cies, where regardless of how shallow is the granular bed 



FIG. 2: Dimensionless rise times for intruders with pr > I for 
all depths, normalized with the rise times of the fastest (dens- 
est) sphere, at each depth. The line is the best fit obtained 
with apr, giving a = 9.84 and b = —1.09. Inset: dimension- 
less rise times for intruders with pr < 1 for different depths, 
normalized with the rise time of the fastest (lightest) sphere, 
at each depth. 



on top of the spheres, they can not segregate if they are 
at the bottom. Indeed, on one hand, at low frequencies 
spheres with pr < 1 can not ascend from the bottom as 
seen in Fign(but spheres with Pr > 1 can). On the other 
hand, at high frequencies spheres with pr < 1 can ascend 
from the bottom (but spheres with pr > 1 can not) . Dur- 
ing this high-frequency experiment tracer particles were 
put at the base of the container, alone or together with 
the light spheres, and the latter emerge but the tracer 
do not (i.e. there is no convection). To explain this ef- 
fect, instead of " void- filling" , we prefer to use a term 
used in fluids and mentioned already by some authors in 
the granular field [2l|: buoyancy. Is this granular effect 
caused by fluidization the missing term for explaining the 
results reported by the group of Chicago [Tsj l and ours in 
Fig^ The answer is no. 

The condition for buoyancy to happen, when neither 
convection nor inertia are present, is pure fluidization. 
Here, we would like to point out that fluidization with no 
convection is achieved if along the granular column there 
isn't a "temperature" gradient; were the particles of the 
granulate move only around their equilibrium positions, 
whether or not the bed crystallizes 

Hence, at high frequencies the granulate fluidizes and 
buoyancy takes place, segregating to the top any light 
intruder buried inside it. Furthermore, at these fluidized 
conditions, where light particles are buoyant, heavy in- 
truders, as in common fluids, must sink. In fig|31we also 
depict sink times as a function of pr for pr > 1. The sink 
curve shows that the heavier is the intruder, the faster it 
sinks (we plot the time it takes for each sphere to sink its 
own diameter). The above concepts explains the RBNP 
predicted by Hong et al 0, later observed by 

Breu et alj25i]. It also explains why Canul-Chay et al |2l| 
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could not observe it in their own experiments |29j . 

To understand the stiU unexplained left part of the 
Mobius et al curve and ours in Fig^ let us explore 
carefully the "convection connection", postulated by the 
group of Chicago ^. We can plot the same data already 
plotted in Fig^ in a different fashion, see the inset of 
the figure. Normalized rise times are plotted as function 
of depth for all relative densities less than one. These 
data were normalized using rise times of tracer parti- 
cles, measured at the same excitation conditions for each 
depth, while no intruders were in the bed. The hori- 
zontal line corresponds precisely to the tracer particles 
rise times normalized to themselves. We note that most 
of the curves are above the tracer line, indicating that 
convection is faster than the ascension of these spheres. 
We suggest then that it is pure convection (neither iner- 
tia because particles are light nor buoyancy because the 
bed is not fluidized at these low frequencies) the mecha- 
nism needed to explain the drop of rise time curves at 
Pr w 1. Looking again to the inset of Fig[51 we see 
that the normalized rise times reasonably collapse into 
the same curve. The fact that they collapse does not 
mean that the spheres rise at the tracer times, as one can 
clearly see in the inset of FigQ] Obviously, although con- 
vection is behind the ascension, their inertia is still non- 
negligible. Therefore, the heavier they are, the harder to 
be dragged to the top by the streaming flux or convective 
cell. This can be understood if we look at the transfer 
of momentum; where both, the "streaming" and "coUi- 
sional" modes participate in the transport of intruders 
having not-negligible masses [28|. 
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FIG. 3: Dimensionless rise (sink) times as a function of rel- 
ative density for only one depth (5 cm), at 50 Hz. In the left 
part of the figure {pr < 1) we note that intruders ascend faster 
the lighter they are, while in the right part {pr > 1) we see 
that the denser sink faster. Inset: rise times as a function of 
relative diameters for pr > 1 (where the solid line represents 
a parabolic fit) and pr < 1 (where the solid line is only to 
help the eye). Both cases are at constant mass (7.1 g and 2.2 
g, respectively), 5 Hz, F = 3, and 5 cm depth. 



Our results plotted in FigQlseem to be in contradiction 
with prior experimental studies conducted by Shinbrot 
et al and Yan et al |2l|. On one hand, Shinbrot et al 
measure convection periods where divergence is seen for 
light intruders. On the other, Yan et al observe that the 
lighter the intruders the faster they sink. In both cases, 
the differences with our experiments lie on the following 
two issues: a) in our experiments we always positioned 
the intruders at specific depths and measure risin g tim es, 
instead of convection periods 0] or sinking times |2lj| , b) 
the size of the bed particles. We are convinced that if 
Shinbrot et al had positioned their intruders at specific 
depths in a granular bed with larger particles, they would 
have observed the non-monotonic rise dynamics we ob- 
serve in FigH Furthermore, Yan et al suggested that in 
order to understand the non-monotonic segregation dy- 
namics shown in Fig^ we need to understand first the 
role of the air pressure gradient acting on the granular 
bed |23|. According to our results, we conclude that this 
suggestion is not correct. Air is indeed an important 
ingredient only in granulates with very small particles 
where reverse or negative |2lj buoyancy is observed. 

Finally, let us remind the reader that all theoretical 
and experimental studies published so far on granular 
segregation agree on the following: the larger the in- 
truder the faster it rises through a vibrating bed. More- 
over, some geometrical approaches conclude that there 
is a critical diameter below which intruders do not as- 
cend 0, Il3j | . The above conclusions are only true when 
the density remains constant as one changes the size 
of the intruders (see for instance, the recent work re- 
ported by the group of Chicago [13 )• Under these cir- 
cumstances, the mass varies with the cube of the diam- 
eter and therefore, due to their greater value, larger in- 
truders ascend faster. However, the opposite behavior is 
observed when the mass (not the density) remains con- 
stant as one changes the size of the intruders, see the 
inset of FiglSJ. In this case, the faster intruders are 
the smaller (in contradiction with the mentioned size 
paradigm). Our above model, that predicts the hyper- 
bolic dynamics for intruders denser than the granulate 
(at constant intruder volume, see Fig|2Jl, predicts also 
the parabolic behavior found for size segregation (see the 
inset of FigE)) at constant mass with p,. > 1 . In the 
inset of FigO we show rising times versus diameter for 
Pr < I (constant mass), verifying in this case what granu- 
lar scientists have concluded in the past: larger intruders 
rise faster. 

We report experimental results that shed definitive 
light to explain the fundamental aspects of the fascinat- 
ing BNP. We conclude that there are only three physical 
mechanisms behind the segregation of a large or small, 
heavy or light, particle in a vibrated granulate: inertia, 
convection and buoyancy (sinkage). The first two are al- 
ways present at F > I and low frequencies, where a non- 
monotonic ascension dynamics is observed as one changes 
the relative density of the intruders. Inertia dominates 
when Pr > 1, and convection does it when p,. < 1. Segre- 
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gation, by buoyancy or sinkage, is present at excitation 
conditions where the granulate is fluidized with no con- 
vection (r > 1, small amplitudes and high frequencies). 
Finally, when intruders have the same mass but diam- 
eters are varied keeping their relative densities greater 
(less) than one, the smaller (larger) they are, the faster 
they segregate to the free surface of the bed. 
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